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We create an ultracold-atoms-based cavity optomechanical system in which as many as six dis- 
tinguishable mechanical oscillators are prepared, and optically detected, near their ground states 
of motion. We demonstrate that the motional state of one oscillator can be selectively addressed 
while preserving neighboring oscillators near their ground states to better than 95% per excitation 
quantum. We also show that our system offers nanometer-scale spatial resolution of each mechanical 
element via optomechanical imaging. This technique enables in-situ, parallel sensing of potential 
landscapes, a capability relevant to active research areas of atomic physics and force- field detection 
in optomechanics. 



Cavity-enhanced interactions between light and the 
motion of a mechanical element form a new structure 
for quantum-sensitive displacement measurements [TJ |2] , 
the generation of non-classical light ^ , the storage [H |5] , 
control [6l [7] and transduction [51 19] of quantum infor- 
mation, and the optical cooling of an oscillator to its 
motional ground state [TOl|TJ. To date, research efforts 
have focused almost exclusively on coupling a single me- 
chanical element to light. 

Cavity optomechanics with multiple moveable ele- 
ments promises new possibilities, including optomechan- 
ical networks for quantum information processing [9l [121 
[T3] , the entanglement of macroscopic objects [HI [15], 
and experimental studies of light-mediated oscillator- 
oscillator interactions [16] . Recently, researchers have be- 
gun exploring multi-oscillator optomechanics, employing 
pairs of nanomechanical oscillators p!7H2Q] . These stud- 
ies have explored coupling of the mechanical modes via a 
shared evanescent optical field or, in the case of Ref. [20 , 
a shared microwave- frequency resonator. Certain clas- 
sical effects were demonstrated, including the synchro- 
nization of motion, and the hybridization of modes into 
optically bright and dark states. 

In this work, we extend these results by observing 
two distinctly quantum-mechanical features of an atoms- 
based realization of cavity optomechanics, in which sev- 
eral distinguishable mechanical elements, prepared near 
their motional ground states, are arrayed within the same 
optical resonator. First, we demonstrate that the near- 
ground-state motion of as many as six oscillators can be 
simultaneously measured with quantum-limited sensitiv- 
ity. Second, we show that the motional state of a chosen 
oscillator can be changed without disturbing the near- 
quantum motion of neighboring oscillators. 

Previous experiments investigating cavity optome- 
chanics through the motion of atomic ensembles have 
used gases trapped within a single-color standing wave 
of light, formed by driving a Fabry-Perot cavity on one 
of its TEMoo resonances [2TH23| . In such a potential, 
atoms at every lattice site oscillate at the same mechan- 



ical frequency. This mechanical degeneracy allows the 
motion of the entire array of trapped ensembles to be 
treated as that of a single mechanical element. 

Here, instead, two optical trap beams of distinct colors, 
labeled A and B, each resonant with a separate TEMqo 
mode, are injected into a Fabry-Perot cavity. Neglect- 
ing the variation of the cavity-field intensity in the di- 
rections transverse to the cavity axis, the resulting one- 
dimensional optical potential is of the form 

U = UaX sin^(/cA^) + % X sin^(/cB^), (1) 

where U and k are the optical potential trap depth and 
wavenumber, respectively, and z refers to the spatial lo- 
cation. In this superlattice, the potential curvature and, 
hence, the mechanical oscillation frequency of trapped 
atoms vary among neighboring lattice sites. Upon dis- 
tributing atoms among several adjacent lattice sites, the 
collective center-of-mass vibrations of the population at 
each occupied site represent distinct mechanical elements 
that can interact optomechanically with a cavity-probe 
field. 

In our experiment, a microfabricated atom chip pro- 
duces an ultracold gas of ^^Rb atoms that is magneti- 
cally trapped within the mode volume of a Fabry-Perot 
optical cavity of half-linewidth = 27r x 1.82 MHz [22] . 
The position of the atomic gas along the cavity axis, 
Za, is controlled by varying the magnetic trap's param- 
eters. Tunable radio- frequency-induced evaporation al- 
lows us to reduce the atom number and also the spatial 
extent of the magnetically trapped atoms by lowering 
their temperature. The remaining atoms are then trans- 
ferred from the magnetic trap to the optical superlattice 
formed by lattice beams A (862 nm) and B (843 nm) (Fig- 
ure [T]). The lattice depths /7a, b are chosen so that the 
mechanical resonance frequencies for each of the single- 
color standing waves that form the superlattice are given 
as a;A,B = \/2^A^t^A^M = 27r x (127kHz, 128kHz), 
where m is the atomic mass (Supplementary Informa- 
tion, Figure S.l). 

Information about the resulting array of oscillators is 



2 



Probe 



843 nm 



862 nm 



IaIa'a", 




Mirror 



s 

. Spatial location ^ , 



Mirror 




f, f. 

Oscillation Frequency 

FIG. 1. Experiment schematic: Superlattice potential. Three 
light beams are injected into an optical cavity (mirrors - grey): 
the probe (dashed blue), and lattice beams A (862 nm - solid 
red) and B (843 nm - solid purple). Atoms inserted into the 
cavity are trapped in potential minima of the superlattice 
formed by spatial beats between each lattice beam's standing- 
wave optical potential (red-purple gradient). The dissimilar 
wavelengths of lattice beams A and B causes atoms at neigh- 
boring lattice sites to resonate at distinct frequencies, which 
is detected via optomechanical interactions with the probe's 
standing- wave intensity (intracavity dashed blue). Since each 
potential minimum overlaps with a different probe intensity 
gradient, the observed response varies in amplitude among 
neighboring lattice sites. 



then extracted via optomechanical interactions with a 
third light beam, the probe beam, resonant with a sep- 
arate TEMoo cavity mode and detuned Aca = —40 GHz 
from the atomic D2 transition (780 nm). At this detun- 
ing, atoms act primarily as a refractive medium with a 
single-atom cavity-QED coupling rate of = 27r x 13.1 
MHz. The displacement, Zi, of each collective atomic os- 
cillator from its trap center, z^, consequently modulates 
the cavity resonance frequency, cjc, and hence the probe's 
phase, by an amount 



(2) 



where kp is the probe wavenumber, Ni is the number 
of atoms in each collective oscillator, and sm{2kp Zi) in- 
dicates the local probe-intensity gradient ^2 E] • The 
cavity's optical emission spectrum thus provides the mo- 
tional spectrum of the optically sensed oscillators. Probe 
light transmitted through the cavity is detected using a 
balanced heterodyne receiver with an overall efficiency 
of 0.13 ± 0.02 for detecting cavity photons. Atoms are 
then released from the optical trap and the cavity probe 
light is again measured. This second measurement is sub- 
tracted from the first to identify and eliminate technical 
noise. To avoid excessive probe-induced heating of the 
atomic gas, data recordings are limited to 10 ms. 



We begin by characterizing the distribution of po- 
tential minima of the optical superlattice. In a single 
experimental cycle, a small cloud of ~ 1000 atoms is 
loaded into the superlattice at an initial position Za cho- 
sen within a range spanning a full superlattice period 
7r/(/cA — ^b) = 19.4 /im. The cavity is then probed at 
a light intensity corresponding to an average intracav- 
ity photon number of 2.9. From the recorded optical 
heterodyne signal, we extract the spectrum of motional 
sidebands imparted by the atomic motion, characteriz- 
ing the mechanical structure of the intracavity atomic 
ensemble. 

These phase-modulation spectra, shown in Figure [2^, 
display a set of discrete mechanical resonances gener- 
ated by atoms trapped at different locations within the 
superlattice. Indeed, the relation between Za and the 
observed mechanical resonances provides a detailed spa- 
tial mapping of the superlattice potential minima, one 
which agrees closely with the predicted form of the po- 
tential. The motional sidebands also vary strongly in 
their strength owing to the site-to-site variation in the 
linear optomechanical coupling strength (Eq.[2|. Indeed, 
the probe-light intensity at each superlattice trap site can 
be inferred from the mean shift of the cavity resonance 
frequency by atoms loaded into that site (Figure ^p)- 
As expected, strong motional sidebands are observed for 
atoms placed where this intensity shows large gradients. 

The precision with which the center location of each 
potential minimum can be imaged depends on the spac- 
ing between adjacent lattice sites, a ~ 420 nm, the ini- 
tial atomic ensemble's spatial full-width at half maxi- 
mum (FWHM), (Jw, and the precision with which Za is 
controlled, <Jz^. For infinitely narrow initial atomic dis- 
tributions with arbitrarily precise initial location, atoms 
inserted within ±a/2 of a potential minimum will be 
trapped at the same location, yielding an effective spatial 
resolution of a. For initial atomic ensembles of FWHM 
c^tot ^ CL^ due to either large or randomly fluctuating 
Za^ each oscillator's center location will be measured with 
a precision of atot = y^cr^ + a^^ . For the region consid- 
ered in Figure |2]:3, the spatial resolution limit is 420 nm. 

Having characterized the optical superlattice, we con- 
struct an array of distinguishable oscillators by inserting 
a spatially broad atomic ensemble (atot ^ 3.1 /im) at 
location 11.7/im (see Figure [2|i,c), where the frequency 
difference between neighboring mechanical resonances 
exceeds the mechanical resonance frequency widths, F 
(FWHM). Atoms are transfered to the many lattice sites 
overlapping with the initial ensemble, resulting in several 
distinguishable mechanical elements composed of 400 to 
800 atoms each. We record the distribution of motional 
sidebands imprinted on the cavity-resonant probe photon 
spectrum, n{uj) (Supplemental Information, section S.l). 
As many as six distinct oscillators are observed (Figure 
|3|. In fact, for these data, atoms occupy at least eight 
neighboring sites of the superlattice, but two sites are not 
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FIG. 2. Distribution of mechanical resonances along the superlattice. (a) Phase-modulation power density as a function of 
frequency and loading position of the initial atomic sample. Discrete mechanical responses appear as bright pixels. They 
closely follow the transparent gray curve, which interpolates between the expected discrete mechanical resonances. Residual 
technical noise appears as horizontal bands in the figure {e.g. near 40 kHz), (b) Measured shifts (purple dots) in the probe's 
cavity resonance are due to dispersive probe-atom coupling at each recorded superlattice position, overlaid with the modeled 
distribution of cavity frequency shifts (red solid line) . Regions of maximum gradient are correlated with the brightest responses 
in (a), (c) Mechanical responses measured with finer steps between loading positions (within dashed box in (a)). Measured 
spectral responses across a selected frequency and position (transparent gray lines) are shown adjacent to the central figure (data 
- purple dots; fits - red solid lines). Fits yield a spatial and spectral resolution of atot = 0.42 zb 0.02 /xm and F = 27r x (1.2 it 0.1) 
kHz, respectively, for the oscillator at frequency cj = 27r x 114.5 kHz. 



detected due to minimal linear optomechanical coupling 
at their locations. 

We demonstrate the quantum nature of these oscilla- 
tors using a parallel measurement of their Stokes emis- 
sion asymmetry. Indeed, the asymmetry between the 
rate of red-detuned (Stokes), n'f\ and blue-detuned 
(anti-Stokes), nf'\ scattered photons, near the mechan- 
ical frequency of element i, reveals that mechanical ele- 
ment to be near its quantum mechanical ground state. 
That is, when the cavity is probed on resonance, this 
asymmetry quantifies the mean phonon occupation to be 
= nf^/[fi[^ — nf^). Mean energy inputs from the 
atoms' average thermal occupation {pf^^) and from op- 
tomechanical interactions with the probe {pf^^^)-, in equi- 
librium with mean energy decay via mechanical damping, 
sum to set Ui'. = z^|*^^ + . Measured occupation 
numbers in our six-oscillator array range from 0.9 to 2.2 
(Figure [3|. These occupation numbers include the per- 
turbative action of the probe and hence constitute upper 
bounds on each oscillator's native motional state, 
As highlighted by the inset of Figure [3) by varying the 
spatial extent of the initial atomic ensemble, the number 
of distinct mechanical elements can be controllably tuned 
from one to six. 



In addition, we demonstrate that one targeted oscilla- 
tor can be driven coherently while preserving the near- 
quantum motion of neighboring oscillators to better than 
95% per added quantum. Oscillators a and (3 (Figure [3| 
are prepared with ~ 750 atoms each. We apply a nearly 
common force modulation to both oscillators by tem- 
porarily modulating /7b, but we tune the modulation fre- 
quency to coincide with one targeted oscillator's mechan- 
ical resonance frequency. Each oscillator's squared dis- 
placement, and hence its mechanical energy, is mapped 
onto the probe's phase-modulation power density, which 
we record at the cavity output (Supplemental Informa- 
tion, section S.2). 

Figure |4] shows the measured simultaneous mean 
phonon occupation of both oscillators. The plotted occu- 
pations include energy provided by the drive in addition 
to the thermal and probe-induced phonon occupation. 
Linear fits to both series of drives indicate a base occu- 
pation of — 1-3. Moreover, for each quantum 
of collective motion added to oscillator a or /3, fewer than 
0.05 and 0.04 phonons are added to the other oscillator, 
respectively. Each oscillator is therefore isolated from its 
neighbor by at least 95% per added mechanical quantum. 
Our findings can be reframed in the context of quantum 
information processing: each oscillator can mechanically 
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FIG. 3. Stokes (dashed red) and anti-Stokes (solid blue) side- 
band spectrum of light exiting the cavity as a function of 
detuning from the probe carrier frequency. Atoms are loaded 
into eight neighboring lattice sites, labeled sequentially by 
position (Greek letters). Only six motional sidebands pairs 
are visible because two oscillators (S and e) are positioned 
where the linear optomechanical response is minimal. Phonon 
occupations are obtained from the integrated Stokes and 
anti-Stokes fluorescence in narrow frequency bands around 
each mechanical resonance (light green): z^c^ = 1.8 =b 0.2, 
ly^ = 2.2 db 0.2, = 1.0 ± 0.1, = 0.9 db 0.1, z/^ = 1.3 db 0.2, 
and 1^0 — 1.8 =b 0.4. The narrow response at 99 kHz is an op- 
tomechanically amplified technical noise spike on the probe. 
(Inset) Demonstration of experimental control over the num- 
ber of constructed mechanical elements. Initial atomic en- 
sembles of narrow (green, top), medium (orange, middle) and 
wide (purple, bottom) spatial extent are prepared to produce 
arrays of 1, 3 and 6 oscillators, respectively. Traces are offset 
for clarity. 



store information quanta for a time l/T ^ 100/is, infor- 
mation that can be retrieved via optomechanical interac- 
tions. Multi-oscillator optomechanical systems thus hold 
potential as quantum registers. 

An important application of cavity optomechanics is 
precise force sensing ^5]. Having an array of mechan- 
ically distinct oscillators allows one to make such force 
measurements simultaneously at many spatial locations, 
using frequency multiplexing to read out the array of 
force sensors with a single cavity-optical output. This can 
serve as a means to measure minute and sharply varying 
force fields, such as short-range gravitational forces and 
Casimir-Polder forces. Additionally, the use of optome- 
chanics to characterize our superlattice potential (Fig- 
ure [2|, a method which could be termed "mechanical 
resonance imaging," can be broadly applied to directly 
measure any potential landscape. Sub-wavelength-scale 
potential variations could be resolved given our system's 
nanometer-scale control over the initial atomic cloud's 
size and location. This technique has direct relevance 
to the many-body physics of quantum interactions [26J, 
Anderson localization ^271 [28] , and studies of particle dy- 



FIG. 4. Selective addressing of a single intracavity oscillator. 
Energy is deposited into oscillator a (green dots) and /3 (blue 
squares), respectively, by coherently modulating the intensity 
of lattice beam B. The points correspond to the measured 
phonon occupation under different modulation strengths of 
Ub- The lines are fits to the data. From the slopes and their 
uncertainties, the phonon occupation of the undriven site in- 
creases by less than 0.04 (blue) and 0.05 (green), within a 97% 
confidence interval, for every phonon added to the driven site. 
(Inset) Measured phase-modulation power density for two 
distinct drives at (purple dashed, orange solid). Arrows 
point to the corresponding data points. 



namics in different optical superlattice structures p9] . 
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SUPPLEMENTAL INFORMATION 

S.l PHONON OCCUPATION - PROBE PHOTON 
SPECTRUM 

For the case of an array of mechanically distinct oscil- 
lators probed on cavity resonance, the intracavity probe 
photon spectrum is given by 



n{uj) 



(SI) 



rf(^. + i) 



(c^-a;,)' + (r^/2)2 (a; + a;,)' + (r^/2)2 

where F^ and Ui symbolize the mechanical full-linewidth 
and oscillation frequency of oscillator respectively. The 
term Ci represents the optomechanical cooperativity of 
oscillator i. The cooperativity parametrizes the rate of 
information exchange between the probe and the oscil- 
lator relative to the rate of information decay by opti- 
cal and mechanical damping. For our experiment, C was 
near unity for each oscillator. In the unresolved-sideband 
limit {k, ^ uji)^ the phonon occupation of oscillator i is re- 
lated to the integrated photon number under the Stokes, 
n-^\ and anti-Stokes, n-^\ mechanical sidebands: 



nl = 



^Qr,(^, + i), 

n) — n) ^ 



(S2) 
(S3) 

(S4) 



We extract n(a;) from the power spectral density 
recorded by our heterodyne photodetector, S\^Qt{uo)^ 
based on the relation 



^het(^) = 5'sN [1 + edet(^(^) + no(cj))/2] , 



(S5) 



where 6'sn = ^lo^^o is the mean shot-noise spectrum, 
with Plo the detected local-oscillator optical power and 
hjjQ the probe photon energy, edet = 0.13 ± 0.02 is 
the overall probe cavity photon detection efficiency, and 
no(co') is the spectrum of technical noise present on the 
probe. Since our probe is observed to be nearly shot-noise 
limited at frequencies near the mechanical resonances, we 
make the approximation no(c^) — 0. We integrate n{uj) 
near each observed Stokes and anti-Stokes mechanical 
sidebands separately to obtain fi[^ and nf^ , respectively. 
The phonon occupation of each oscillator i is then calcu- 
lated from Eq. IS4 



S.2 PHONON OCCUPATION - PROBE 
PHASE-MODULATION POWER 



The phonon occupation of oscillator a and /3 under 
applied force modulations is determined from the probe's 
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FIG. SI. Schematic of experimental optical cavity. Lattice beams A (red) and B (purple), and the probe (blue) are independently 
monitored in transmission of the optical cavity. By adjusting their measured transmission to a desired level, [/a, Ub and the 
intracavity probe intensity can be independently controlled. Lattice beam B's transmission signal is also used to feedback and 
lock the cavity's length. A milliwatt of 780-nm light circumvents the cavity and is used as the local oscillator (orange) for 
balanced heterodyne detection of the pico- watt- level probe. Recorded probe information is used to lock the local oscillator path 
length (quadrature lock), to control the probe-cavity detuning (probe lock) and for data analysis. Probe frequency control is 
accomplished by placing weak- intensity frequency markers on the probe and feeding back based on the lower frequency marker. 
Acronyms: AO = acousto-optical modulator, EO = electro-optical modulator, APD = avalanche photo-diode. 



phase spectrum, n^^/^^(co'), which is obtained from the 
recorded phase- modulation power density, Sl^J'^\uj): 



o(^/2) 
^het 



(7r/2) 
SN 



€det 

2 



n 



(7r/2) 



(S6) 



where S^^'^^ = is the mean phase- modulation 

shot-noise spectrum recorded, and nQ^^'^\uj) is the spec- 
trum of phase-modulation technical noise on the intra- 
cavity probe field. As before, we observe the probe spec- 
trum to be shot noise dominated near mechanical side- 
bands, and consequently approximate nQ^^'^\uj) to be 
0. For a cavity-resonant probe, the mechanical motion 
modulates only the probe's phase, and hence both the 
Stokes and anti-Stokes motional sidebands are mapped 
onto n(^/2)(cj): 



K 



CiYl{2vi + 1) 



(S7) 



where F^, a;^, Ci and vi are the mechanical linewidth, 
mechanical resonance frequency, optomechanical cooper- 



ativity and mean phonon occupation of oscillator z, re- 
spectively. 

We experimentally determine parameters relevant for 
establishing vi in the presence of an applied force. The 
thermal phonon occupation of both oscillators, z^^*^^ ~ 



,(th) 



is first obtained from time-of-flight ther- 
mometry of the atomic gas [2T in the absence of a probe 
and an applied drive. Mechanical properties are then 
determined by probing oscillators a and /3 in the ab- 
sence of a drive. Fits of the observed motional sidebands 
in n^'^l^^iyS) are used to extract F^, uji and Ci, where 
i = a, /3. Each oscillator's cooperativity is derived from 
the peak height of its motional sideband: 



ni"/'\t^i) = 4Ci(2Pr' + + 1) 



(th) 



(S8) 



which includes the unresolved-sideband approximation 
(^K ^ uji) and the probe's radiation-pressure shot noise 
backaction on the mean phonon occupation (z^j^^^^ = 

Ci /2)m- 



